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We summarize our recent results on the ground state energy of multi-polaron systems. In 
particular, we discuss stability and existence of the thermodynamic limit, and we discuss 
the absence of binding in the case of large Coulomb repulsion and the corresponding 
binding-unbinding transition. We also consider the Pekar-Tomasevich approximation to 
the ground state energy and we study radial symmetry of the ground state density. 
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1. The Frohlich Hamiltonian for a single polaron 

The large polaron was first considered by H. Frohlich in 1937 as a model of an 
electron interacting with the quantized optical modes of a polar crystal.^ In suitable 
units, its Hamiltonian is 

ijW(a) =p 2 -VoXx)+ / a t (k)a(k)dk (1) 

JUL* 
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with 

1 f 1 

y>(x) = -=- I -r- [a(k) cxp(ik • x) + h.c] dk . 
V2 7T Jr3 \k\ 

This Hamiltonian can be defined as a self-adjoint and lower semi-bounded operator 
in the Hilbert space L 2 (M 3 ) ® J 7 , where T is the bosonic Fock space over L 2 (R 3 ) for 
the longitudinal optical modes of the crystal, with scalar creation and annihilation 
operators a'(k) and a(k) satisfying [a(k), a^(k')] = <5(k — k'). The momentum of an 
electron is p = — iV, and the coupling constant is a > 0. The ground state energy 
of HW(a), 

E {1) {a) =infspe Cj ff (1) (a), 
has been studied in detail and we summarize the following properties. 

(i) For all a one has the uppeJ^- 17 * 18 ! and lowei^ bounds 

—a — a /3 < E (1) { a) < —a . 

As a consequence, E^ (a) ~ ~a for a small. 

(ii) Using a product function Pekap 2 ^ showed that 

E^(a) < -C P a 2 

for all a. Donsker and VaradharP showed that this bound is asymptotically 
correct and Lieb and Thomas^ obtained the error estimate 

E (1 \a) > -C P a 2 - const a 9 ' 5 

for large a. Here, Cp — 0.109 is the number determined by Pekar's variational 
problem for the electron densitypSl 

Cp =mU m»d*-[[ i^)i 2 iwi 2 , xrfy:M2 = 1 \, (2 ) 



The minimizing ip is unique up to translations and a trivial phased 
(iii) There is a representation for E^(a) in terms of path integrals.- In terms of 
the partition function z£\a) = Tr exp ( - TijW(a)), one has E^\a) = 
— limT_i.oo T~ Y log Z^\a). (Strictly speaking, zfjp(ot) does not exist because 
of the translation invariance of H^^a) and the infinite number of phonon 
modes. These technicalities can be handled by inserting appropriate cutoffs, 
to be removed at the end of the calculation! 28 ^ 33 !) After one integrates out the 
phonon variables, Zjp(at) has the functional integral representation 



Z^\a) = J dfiexp ~ J J 



T -"^l dtds 



lo Jo |x(t) - x(s)| 



(3) 



where d\i is Wiener measure on all T-periodic paths x(t). (In physics notation 
d[i = exp(— J Q x(i) 2 dt) (ipath. Strictly speaking, t — s has to be understood 
modulo T, but this is irrelevant as T — ¥ oo.) 
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2. Multi-polaron systems and their thermodynamic stability 

The Hamiltonian for N polarons is 

H ( u N \a) = Ef=i (p? - V*<P<&)) + /». « f (k)a(k) dk + UE 1 < t<J < N ^ (4) 

and we denote its ground state energy by 

E ( ^\a) = infspecfl£ W) (a). 

This operator acts in the Hilbert space L 2 (R 3JV ) ® F. We ignore Fermi statistics 
for the electrons, because its imposition changes things only quantitatively, not 
qualitatively. The Coulomb repulsion parameter U > is equal to e 2 . Frdhlich's 
derivation^ of H^\a) implies that U > 2a, and this is crucial for thermodynamic 
stability, as we shall see. 

We consider the question of the existence of the thermodynamic limit for a multi- 
polaron system in the ground state. For large N, physical intuition suggests that 
E^\a) r~j — const N. This supposition is known to be false if U < 2a. Indeedp' 
even with the Pauli principle, E^f\a) ~ — const N 7 / 3 when U < 2a. Absent the 
Pauli principle, E^\a) would behave even worse, as — const TV 3 . It is also knowifEQ 
that E\f\a) > — const N 2 if U > 2a. The latter bound ought to be — const iV 
instead, and this is indeed the statement of the following theoremP^ 

Theorem 2.1 (Thermodynamic stability for U > 2a). For any U > 2a > 

0, there is a constant C(U,a) such that for all N > 2, 

E { uHa) > -C(U,a) N. 

Our lower bound on N~ 1 E < jj f \a) goes to — oo as U \ 2a, but we are not 
claiming that this reflects the true state of affairs. Whether limjv->oo N~ 1 E^ (a) 
is finite or not remains an open problem. There are results on this question in the 
Pekar-Tomasevich approximation!^-!] 

The linear lower bound from Theorem 12. 1\ together with the sub-additivity of 
the energy EMI [ e ; 

E ( u N+M \a)<E^\a) + E^\a), (5) 

implies: 

Corollary 2.1 (Thermodynamic limit for U > 2a). For any U > 2a > 0, 

limAT^oo N^ 1 E^j'\a) exists. 

3. Binding and non-binding of multi-polaron systems 

The binding of polarons, or its absence, is an old and subtle topic. For some time 
the bipolaron binding energy AEu(a) — 2E ( - 1 \a) — E^\a) was thought to be zero 
for all U > 2a, on the basis of an inadequate variational calculation, but it is now 
knowrPto be positive for some U > 2a. The first question we address is whether 
AEjj(a) = for U sufficiently large. It is understood that the effective interaction 
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induced by the phonon field for two polarons at large distances d is approximately 
Coulomb-like —2a/d, but this alone does not preclude binding. (The reason for 
2a ■ distance -1 can be seen from the TV-polaron analogue of ([3]), 

i-T r T __|*-.| j. . pT 



dt 



, ; •• I *i (t ) ~ Xj (t ) | 

where d/z^) is Wiener measure on all T-periodic paths (xi (t), . . . , xjy(t)). There is a 
factor a/2, but the pair (i, j) appears twice, and the integral / R e~' t_s ' ds = 2.) The 
known existence of bipolarons for some U > 2a is an effect of correlations. It is a pri- 
ori conceivable that correlations lead to an effective attraction that is stronger than 
Coulomb at large distances. If it were, for example, equal to (2a/ d) log(log(log(d))), 
then this minuscule perturbation of Coulomb's law, which would be virtually unde- 
tectable by a variational calculation, would result in binding for all U. The absence 
of binding is a problem that has resisted a definitive resolution for many years. We 
proved 

Theorem 3.1 (Absence of binding for TV polarons). For any a > there is 
a U c (a) < oo such that for all U > U c (a) and all N > 2 one has 

E^ ] (a)^NE^(a). (6) 

In particular, for TV — 2 we show that E^'(a) = 2E^(a) provided U > 2Ca 
with C — 26.6. The constant 26.6 vastly exceeds the current, non-rigorous estimates 
of about 1.15,^2121] so it is an open problem to find a more accurate rigorous bound. 
Somewhat better bounds are known in the Pekar-Tomasevich approximation!^ 

While our bound for U c (a) is linear in a for large a, we have not achieved this 
linear bound for small a and this remains an open problem. 

Theorem 13 . 1 1 says that 



C/ c (JV) (a) = inf |Z7 > : E^^a) = NE^(a) for all U' > If} 

is finite and bounded uniformly in TV. For any U > ui N ^ (a) and any state ^ 
(tf|f#*°(a)|tf) >NE^(a)\\^f+(u-U^(a)^ (*| ^ Ix^-Xjl" 1 ^ . (7) 

i<j 

This is a quantitative estimate of the energy penalty needed to bring two or more 
particles within a finite distance of each other. In particular, it implies that for U > 
Uc (oi) there cannot be a normalizable ground state, even in a fixed momentum 
sector. Inequality ([7]) is not only true for our bound on U c (a), but also for the 
(unknown) exact value of Uc (a). 

For U in the range 2a < U < U c (a), there are bound states of an undetermined 
nature. Does the system become a gas of bipolarons, or does it coalesce into a true 
TV-particle bound state? If the latter, does this state exhibit a periodic structure, 
thereby forming a super-crystal on top of the underlying lattice of atoms? This is 
perhaps the physically most interesting open problem. While particle statistics does 
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not play any role for our main results, the answer to this question will crucially 
depend on particle statistics (Bose or Fermi) 



4. Binding unbinding transition 

We now discuss the behavior of the iV-polaron radius as the repulsion parameter 
U approaches from within the binding regime a critical value where ./V-polarons 
cease to be bound. Does the iV-polaron radius in this limit increase towards infinity 
or does it remain finite? Verbist, Peeters and Devreese^ proposed a 'first-order' 
transition, that is, the Coulomb repulsion jumps discontinuously from a positive 
value to zero and the radius, too, jumps discontinuously. We prove this rigorously 8 
under the assumption that the critical value is strictly bigger than 2a. This is known 
to be satisfied for large aP3 

Besides the ground state energy E^ (a) of the Hamiltonian (a) we also 
need the minimum break-up energy 



(a)= min (e^ (a) + E$ ~ n) (a)) . (8) 

Kn<N—l \ U U J 



Note that it is always the case that E\j '(a) < E\j (a). 

Theorem 4.1 (Upper bound on the 7V-polaron radius). For any N > 2 and 

e > there is a constant C e (N) > such that for all < 2a(l + e) < U with 
E\f\a) < E^\a) and all states * 



, U-2a(l + e) (^\E^(a)-H^(a) 
* > > 



C e (N)(l + U/a) E^\a)-E^ N \a) 



(9) 



Since (JSJ holds for all "J, it can be reformulated as an operator inequality The 
bound is non-trivial only for states * with (^\H^\a)\^!) < Eff' (a)\\^\\ 2 , how- 
ever, which exist since Eff {pi) < E^\a) by assumption. For approximate ground 
states, that is, states satisfying 

(¥|ffg%)|*) < (1 - 0)E^ N \a)m 2 + 9E^ N \a)\\n 2 
for some 9 > 0, (j9]) gives the uniform lower bound 

1 



U-2a(l + e) „ , 
' - C e (N)(l + U/a) 11 1 



2 



LdXj^j |Xj X 

which means a uniform upper bound on the radius of the multipolaron system. This 
bound depends only on the value of 9 and does not explode as U approaches the 
critical unbinding value. 

A similar phenomenon was shown by T. and M. Hoffmann-Ostenhof and B. Si- 
morP^l for a two-electron atom and the proof strategy of Theorem 14.11 applies to 
that problem as wellP The overall lesson is that this kind of discontinuous binding 
will occur whenever the net repulsion at large distances falls of slower than r~ 2 . 
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5. The Pekar Tomasevich approximation 

The Pekar-Tomasevich approximation to the ground state energy (a) consists 
in minimizing fif\Hy (ot)\^) only over "f's of the form ip(&&, where ip £ L 2 (M 3Af ), 
$ G J-, and both ip and $ are normalized. For N — 1 we obtain functional ((2]). In 
the iV-polaron case this approximation leads to the minimization of the following 
Pekar-Tomasevich functional for normalized functions ip on M. 3N , 

N 



' N Xj-X,- ././k3vB3 x-y 



. 1 Ju 3N — ' 
i=i K %<] 



where dX = rTfc=i ^ x &- The density p^p of ip is defined as usual by 



A' 



p^,(x) =y~] |^(xi, . . . ,x, . . . ,xat)| 2 dxi • • -dxi ■ ■ -dx N 

~{ JR3(«-1) 

with x at the i-th position, and dx^ meaning that dxj has to be omitted in the prod- 
uct n^—j rfxfc. Since the Pekar-Tomasevich functional is the result of a variational 
calculation, its energy gives an upper bound to the ground state energy E^(a). 

The Pekar-Tomasevich minimization problem has been studied in great de- 
tail! 1|2 | 7| E 9 | 11 | 2 ° | 21 1 2 -^ In particular, the analogous statements of our Theorems 12. 1[ 
13.11 and 14.11 remain valid for this model. 

Here, for the sake of simplicity, we restrict ourselves to the case of a bipolaron, 
so that the energy functional becomes 

£ u[ip] = JJ R3xR3 (|V X ^| 2 + |V y ^| 2 + ^H 2 ) dxdy 

, rr P^.(x)p^(y) , , 

• JJr 3 xR 3 F^yl ax ay . 



-a 



The electron coordinates are x and y and the electron spin does not appear explic- 
itly, except that ip is symmetric for the ground state, which is a singlet state. 

There is a considerable literature on the subject of rotation invariance of the 
bipolaron energy minimizer, usually formulated in the language of 'one-center bipo- 
laron versus two-center bipolaron' EltlfiM] T/ ne analyses are all based on variational 
calculations. While there seems to be general agreement that the one-center bipo- 
laron has the lower energy, it is not completely clear that a more sophisticated 
variational treatment will preserve rotational symmetry, especially near the value 
of U where the bipolaron ceases to be bound. 

The value of U determined by physical electrostatic considerations is always 
U > 2a. Nevertheless, one can consider the mathematical question for small, but 
positive U and ask whether there is a possible lack of rotational invariance in that 
case. After all, a rotating object like the earth becomes oblate even for the smallest 
amount of rotation. 

Our theoreirP says that for the bipolaron rotational symmetry is not broken for 
small U: 



Theorem 5.1. There is a v s > such that for all U < v s a the minimizer of £u 
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is unique up to translations and multiplication by a constant phase. In particular, 
after a translation it is rotation invariant, that is, ip(TZx,TZy) = ip(x., y) for any 
x, y e K 3 and any K € 0(3). 

Our proof of Theorem 15.11 is perturbative in nature. It uses crucially a non- 
degeneracy statement about the single polaron functional ((2]) due to LenzmannP^ 

We also show that for certain values of U /a the bipolaron equation has a positive 
solution which is not a minimizer. 

It remains an open problem to decide whether the ground state ceases to be 
rotation invariant for U close to the critical value where bipolarons cease to be 
bound in the Pekar-Tomasevich approximation. 
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